cpi uiaii 



Zero-temperature spin-glass freezing in self-organized arrays of 
Co nanoparticles 



^ ! R. Lopez-Ruiz 1 , F. Luis* 1 , J. Sese 2 , J. Bartolome 1 , C. Deranlot 3 and F. Petroff 3 

1 Instituto de Ciencia de Materiales de Aragon - CSIC-Universidad de Zaragoza, 50009 Zaragoza, Spain, and 
. Departamento de Fisica de la Materia Condensada, Universidad de Zaragoza, 50009 Zaragoza, Spain 

2 Instituto de Nanociencia de Aragon, Universidad de Zaragoza, and Departamento de Fisica de la Materia Conden- 
, sada, Universidad de Zaragoza, 50009 Zaragoza, Spain 

C*~) ' 3 Unite Mixte de Physique CNRS/Thales - Route Departementale 128, 91767 Palaiseau Cedex, France, and Universite 
£N , Paris-Sud - 91405 Or say Cedex, France 



G 

O 
O 



> 

CN 

OS 

o 



X 



PACS 75. 50. It - Fine-particle systems; nanocrystalline materials 

PACS 75 . 40 . Gb - Dynamic properties 

PACS 75 . 50 . Lk - Spin glasses and other random magnets 

Abstract. - We study, by means of magnetic susceptibility and magnetic aging experiments, 
the nature of the glassy magnetic dynamics in arrays of Co nanoparticles, self-organized in N 
layers from N — 1 (two-dimensional limit) up to N = 20 (three-dimensional limit). We find no 
qualitative differences between the magnetic responses measured in these two limits, in spite of 
the fact that no spin-glass phase is expected above T = in two dimensions. More specifically, 
all the phenomena (critical slowing down, flattening of the field-cooled magnetization below the 
blocking temperature and the magnetic memory induced by aging) that are usually associated 
with this phase look qualitatively the same for two-dimensional and three-dimensional arrays. The 
activated scaling law that is typical of systems undergoing a phase transition at zero temperature 
accounts well for the critical slowing down of the dc and ac susceptibilities of all samples. Our data 
show also that dynamical magnetic correlations achieved by aging a nanoparticle array below its 
superparamagnetic blocking temperature extend mainly to nearest neighbors. Our experiments 
suggest that the glassy magnetic dynamics of these nanoparticle arrays is associated with a zero- 
temperature spin-glass transition. 



Introduction. — Dense arrays of magnetic nanopar- 
ticles contain the physical ingredients of spin-glasses pQ. 
Disorder in the positions and orientations of the particles 
leads to disorder and frustration of the dipolar interac- 
tions, usually dominant, between their magnetic moments. 
In contrast with "canonical" spin-glasses, the slow mag- 
netic relaxation introduced by interactions coexists with 
the slow magnetization reversal associated with the high 
anisotropy energy barriers. Many experiments performed 
on dense nanoparticulate materials show phenomena, such 
as magnetic aging [2] and the slowing down of the ac 
susceptibility [3] , which are typical of spin-glasses [4HTT] . 
However, some of these phenomena are not exclusive of the 
spin-glass phase The question is, then, whether 

real nanoparticulate materials show a true (super)spin- 
glass phase. 



Experimental studies are often hindered by the lack of 
control over the sample parameters that determine the 
nature and strength of dipolar interactions, such as inter- 
particle distances, spatial organization, etc. This usually 
makes it difficult to know a priori if a particular system 
is expected to show a spin-glass phase. Perhaps the most 
clear-cut situation to discuss the existence of a phase tran- 
sition and its experimental manifestations is offered by the 
study of a single layer of nanoparticles. In contrast with 
three-dimensional systems [2], it is generally accepted 
[151116] that the transition temperature T g vanishes in two- 
dimensions. Results of tempered Monte Carlo simulations 
seem to confirm the same conclusion also for Ising spins 
coupled by dipolar interactions [17) . 

Based on these considerations, our work was aimed to 
elucidate the nature of the glassy magnetic dynamics, i.e. 
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whether it is associated with a superspin-glass phase at 
T g > or if, by contrast, T g = 0, in self-organized 
nanoparticle arrays. For this, we compare results obtained 
on very well-characterized three- and two-dimensional ar- 
rays of Co nanospheres. Previous experiments reveal that 
the superparamagnetic blocking temperature Xb, defined 
as the temperature of the in phase yf ac susceptibility 
cusp, increases as additional layers are deposited on a 
two-dimensional sample [18 1 119 ) . Since the number of 
layers modifies the number of nearest neighbors in the 
nanoparticle array, that result indicates that dipolar inter- 
actions slow down the magnetic relaxation processes. In 
the present study, we have investigated how the number of 
layers modifies the critical slowing down and the magnetic 
aging, properties that are usually associated with the spin 
glass behavior [3D]. Our results show that no qualitative 
changes in these quantities occur as the two-dimensional 
limit is approached. The control over the number of lay- 
ers and their separation has also enabled us to directly 
probe the magnetic correlation length and show that it is 
mainly restricted to a first shell of nearest neighbors and, 
in any case, shorter than what would be expected for a 
conventional spin glass. 

Experimental details. — Samples made of N lay- 
ers of Co nanoparticles with average diameter D ~ 2.6 
nm were prepared by the sequential sputtering of N = 
1,2,3,4,5,7,10,15, and 20 Co and A1 2 3 layers on sili- 
con substrates [H1[T11[2TJ[22]. The particle's shape and 
average size (thus also the average magnetic moment // p 
per particle), as well as the width of the size distribu- 
tion (<td = 0.26-D) are approximately independent of iV 
[19] . Nanoparticles deposited on adjacent layers tend 
to self-organize in a structure that resembles a closed- 
packed hexagonal lattice of nanospheres [22 . The separa- 
tion between the Co layers is determined by the thickness 
^ai 2 3 = 3 nm of the alumina layer. Nearest neighbors 
separations are d nn u ~ 4.6 nm, within a given layer, and 
rfnn,_L — 4.2 nm, between adjacent layers. They correspond 
to dipolar energies E^ip = /J.p/d^ n sa 13 K and 17 K, re- 
spectively. As described in [23] , the anisotropy energy bar- 
rier Uq for the magnetization reversal was estimated from 
ac susceptibility experiments performed under sufficiently 
strong magnetic fields, which dominate over dipolar inter- 
actions. This method gives Uo ~ 430 K. In the same way, 
we estimate an attempt time tq ~ 10~ 13 s, of the same 
order of that found for samples of very small Co nanopar- 
ticles (D ~ 1 nm), prepared by the same technique [18] . 
for which interactions are expected to become neglibible. 
A multilayer with N — 20 layers but a larger interlayer 
separation iAi 2 o 3 = 10 nm, and thus also a much smaller 
interlayer -Edip ~ 1-6 K, was prepared under identical ex- 
perimental conditions. 

Ac susceptibility and magnetization measurements were 
performed with a commercial SQUID magnetometer. 
Samples were rectangular plates with approximate dimen- 
sions 9 x 3 x 0.5 mm 3 . Ac and dc magnetic fields were 



parallel to the plane of the sample to minimize demagne- 
tizing effects. In our study of aging [5J[7], we measured the 
time-dependent relaxation of the zero-field cooled (ZFC) 
magnetization on samples aged, at zero field, for a time t w 
at temperatures T w < T^. In addition, we employed a dif- 
ferent method which consists on measuring magnetization 
curves (zero-field and field cooled (FC), and remanence) 
using the waiting time protocol described in [241125] . 

Results and discussion. — A typical method to 
characterize the spin-glass behavior is by measuring the 
frequency-dependent ac magnetic susceptibility [3-5 r At 
any fixed frequency uj we define a characteristic relaxation 
time t c such that wr c (T) = 1 at T = Tb(ui). For spin 
glasses t c diverges at T g according to a power law, reflect- 
ing the growth of magnetic correlations [4] 

r c = T*|l-T/T g |— (1) 

In Fig. [TJ we plot r c versus the reduced temperature for 
N ranging from 1 to 20 layers. The experimental data are 
compatible with a critical slowing down of the magnetiza- 
tion dynamics at a finite T g . In order to limit the number 
of fitting parameters, we took T g , for each sample, as the 
temperature of the ZFC susceptibility cusp (i.e. equal to 
the 7b corresponding to a typical timescale of the order of 
170 s). The microscopic time scale and the dynamical 
critical exponent zv are found to be nearly the same for 
all samples. The critical exponent is close to typical values 
found for spin glasses [5] as well as for some nanopartic- 
ulate materials [25]. The present results are remarkable 
because it is generally believed that T g = in two dimen- 
sions. Notice however that, as often happens with plots of 
this type obtained for nanoparticles 3,26,27 , experiments 
do not explore the close vicinity of the critical region. For 
this reason, the data are relatively easy to fit. Fits of sim- 
ilar quality can be obtained by scaling all T g 's by a factor 
in between 1 and 0.75. The characteristic t* increases 
then from about 10~ 6 s to 10~ 4 s while, at the same time, 
the exponent zv increases from 7.3 to 14. In fact, if one 
wishes to include also in the analysis the freezing tem- 
perature extracted from ZFC susceptibility data (getting 
closer to T g ), the best fits with Eq. ([TJ are obtained then 
for the largest zv and values (thus also for the lowest 
T g ). Such large zv values are not uncommon in systems 
of magnetic nanoparticles [3 , 27J but they are significantly 
larger than what it is expected for a canonical spin glass 
phase transition (of the order of zv = 7 28J ) . 

An alternative theoretical framework to describe the 
frequency-dependent susceptibility, which seems very ap- 
propriate in the case of a layered material with a markedly 
two-dimensional character, is the activated dynamics char- 
acteristic of glassy systems undergoing a phase transition 
at T g = [HQS]. in the latter situation, the critical slow- 
ing down of t c obeys the following expression 

t c = T cxp (Ejk B T f (2) 
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Fig. 2: (Color online). Log-log plot showing the variation with 
temperature of ln(r c /ro), where r c is a characteristic relaxation 
time extracted from ac susceptibility data and to = 1CP 13 s. 
The lines are fits of the law r c = to exp (E^/k^T)" , character- 
istic of a spin-glass transition at T g = 



Fig. f: (Color online). Critical slowing down of the charac- 
teristic relaxation time extracted from ac susceptibility exper- 
iments. Results are shown for samples with varying number of 
layers N. 



where E a is an effective activation energy and a is a crit- 
ical exponent. As it is shown in Fig. [2] we find a good 
agreement with our data, including also the temperature 
of the ZFC magnetization cusp, for a = 1.3 (to be com- 
pared with a = 3.2 found for 2 — D spin-glasses [TB]) and 
E a gradually increasing with the number of layers from 
345 K up to 471 K. From these frequency-dependent sus- 
ceptibility experiments, we conclude that the nature of 
the slow magnetic dynamics of two-dimensional (i.e. with 
N equal to or close to unity) and three-dimensional (with 
large N) nanoparticle arrays is the same. The description 
based on a zero-temperature phase transition is appeal- 
ing, because it is consistent with the behavior expected 
for a single layer. By themselves, however, these experi- 
ments cannot discriminate between the two alternatives, 
i.e., whether the underlying physics corresponds to the ex- 
istence of a second-order phase transition at a finite T g or 
if, by contrast, T g = 0. 

Aging experiments can shed some light and help decid- 
ing between these two alternatives, since they probe how 
dynamical magnetic correlations grow with time jUIHEI] ■ 
We have carried out two different experiments, which mea- 
sure the magnetic memory effects associated with the 
aging of the sample at a given temperature. In the 
first of these, the quantity of interest is the difference 
AM = M - M w between the magnetizations (ZFC, FC 
or remanent) measured after cooling the sample without 
or with a pause at an intermediate temperature T w < Tb 
[24"] . Results measured for t w = 10 4 s are shown in Fig. [3] 
AMzfc shows a peak centered near T w . If T w is varied, 
the peak shifts accordingly. In addition, the relationship 



[24] AMpc = AM r +AMzpc is fulfilled, showing that they 
are associated with the aging of the sample at T w and not 
with experimental artifacts. 

Figure |3] compares results obtained on a single layer 
N = 1 with those measured on a multilayer made of 
N = 20 layers. The aging was performed at T w = 0.7Tb 
for the two samples. Besides the obvious difference in 
the signal-to-noise ratios, they look qualitatively the same. 
The maximum in AMzfc vs T is just about 25 % larger in 
the case of the multilayer. A first conclusion is, therefore, 
that the magnetic memory induced by aging a nanopar- 
ticle array does not show any abrupt change as the two- 
dimensional (2D) limit is approached. Notice also that, as 
we have seen with the critical slowing down, the analogy 
is not restricted to aging. The FC curves measured on the 
two samples show also the same degree of "flattening" be- 
low Tb, a property that has been considered as a signature 
of the superspin glass phase [25] , 

By gradually changing the number of layers N we can 
study how magnetic correlations grow. In Fig. |4] we show, 
as a function of N, the relative amplitude of the mag- 
netic memory effect AMzfc/Mzfc measured after aging 
at T w = 15.8 K for i w = 10 4 s. Within the droplet picture 
of the spin-glass phase [29] , this quantity is connected with 
the size that domains of correlated spins attain after time 
tw [II]- We see that AMzfc/Mzfc increases rapidly when 
one or two layers are added to a two-dimensional sample, 
nearly saturating as N increases further. The right-hand 
panel of Fig. U shows that, within the relatively large ex- 
perimental uncertainties, AMzfc/Mzfc ^ s approximately 
proportional to the increase in the average number of near- 
est neighbors N± — 6(N—1)/N that is associated with the 
addition of extra layers. The same linear dependence was 
also found for the blocking temperature Q15] . This behav- 
ior suggests that the enhancement in the amplitude of the 
magnetic memory is provided mainly by correlations with 
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Fig. 4: (Color online). Left: Variation with the number N 
of layers of the memory AMzfc /A^zfc measured after aging 
the sample for i w = 10 4 s at T w = 15.8 K (•). Right: Same 
data as a function of the average number of nearest neighbors 
N± = 6(N — 1)/N that a nanoparticle has in adjacent layers 
p] , Lhe lines represent 0.033 + 0.015JVx. 



Fig. 3: (Color online). ZFC, FC and thermoremanence curves 
measured with (red crosses) and without (black solid line) 
pause. The insets show absolute values of AMzfc, AMfc and 
AM T . Left: sample with N = 1; T w was 15.8 K. Right: sample 
with N = 20; T w was 25 K. The waiting time was 10 4 s and 
the applied magnetic field was 10 Oe. 

the first one or two nearest layers. Another result suggest- 
ing that magnetic correlations remain rather short-ranged 
is shown in Fig. [5l There, we compare the magnetic mem- 
ory AMzfc obtained for a single layer of 2.6 nm particles 
with that obtained for a N = 20 multilayer in which the 
interlayer separation is £ai 2 o 3 — 10 nm, i.e. more than 
twice d nn ii ~ 4.6 nm. Within their respective experimen- 
tal uncertainties, these two quantities are found to be the 
same. Also Xb and other quantities agree. It seems then 
that no measurable magnetic correlations are established 
between layers of nanoparticles located 10 nm far form 
each other. 

We also studied the effects of aging using a different 
experimental method, which enables a more quantitative 
determination of magnetic correlation lengths. For this, 
we measured the magnetic relaxation of the ZFC mag- 
netization of a N — 15 multilayer at T w = 15.8 K. The 
sample was first cooled from 100 K to T w in zero field. 
After aging the sample for t w — 10 4 s, a magnetic field 
H was applied and the ensuing magnetization measured 
as a function of time. This method has been applied to 
estimate the number of correlated spins in spin glasses 
[30, 31], and recently applied also to investigate the slow 
dynamics of frozen ferrofluids 26 . Its basic idea is as 
follows. During the waiting time t w , at zero field, mag- 
netic correlations between nanoparticles grow [6]. Typi- 
cal free energy barriers A(i w ) for the flip of iV s (i w ) cor- 



related spins increase also with the age t w of the sys- 
tem. This growth of dynamical correlations reflects itself 
in the appearance of a maximum in the relaxation rate, 
defined as <9MzFc/<91og(i), when the experimental time 
t approaches the age of the system i^f ~ t w (see Fig. 
[5]). A magnetic field H reduces the free energy barriers, 
from its zero-field value to A(t w ) — Ez \H, iV s (i w )], where 
E Z [H, N s (t w )} = n[H, N s (t w )] H and fx is the magnetic 
moment of a "drop" of N s correlated spins. The energy 
shift induced by this Zeeman term effectively reduces the 
" age" of the system according to 

C ff (i7) = C(^.0)ex P {-^l^M} (3) 

therefore shifting the relaxation rate maximum towards 
shorter times with increasing H , as it is indeed observed 
experimentally (Fig. [5]). From a series of experiments 
performed at different fields, ranging from 5 Oe up to 100 
Oe, we have extracted the Zeeman energy Ez [H, iV s (i w )], 
which we plot in the main panel of Fig. [6] The field 
dependence of this energy can be fitted using a quadratic 
function of H, compatible with the following expression 

Ez [H, N s (t w )} = N s (t w )xz FC H 2 (4) 

which was found to agree also with experiments performed 
on frozen ferrofluids Inserting in Eq. (Q} the mea- 

sured FC susceptibility per particle xfc ; we estimate the 
number of correlated spins N s to be approximately 17, i.e., 
rather close to the average number of nearest neighbors 
(12) in a multilayer of nanoparticles [19] , 

Our experimental findings suggest therefore that mag- 
netic correlations achieved after aging the sample at low T 
extend mainly to nearest neighbors. The detailed charac- 
terization of our samples enables us to make a quantitative 
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Fig. 5: (Color online). Magnetic memory for a single layer of 
2.6 nm particles (black thin line) compared with data obtained 
for a N = 20 multilayer in which the interlayer separation is 
^ai 2 o 3 — 10 nm (red thick line). 



comparison of the present results with predictions for the 
growth of correlations in spin-glasses. Theoretical con- 
siderations as well as experiments support the idea that 
correlations grow approximately as a power law of time 

|25|C2ai32|: 



e(t*,T w )/d nn ~(f) 



a(T v ) 



(5) 



where d nn ~ 4.4 nm is the distance to nearest neighbors, 
t* = t w /r(T w ) is a dimensionless timescale, the exponent 
a = 0.17(T W /T g ), and r is the relaxation time of individ- 
ual spins at the given temperature. We have estimated 
t = tq exp {U/ksT) using parameters estimated, as de- 



scribed above, for the noninteracting case: r ~ 10 



-13 



10 4 



and U ~ 430 K [23]. For T w = 15.8 K and t v 
Eq. Ogives 4.2d nn < £ < 7d nn , i.e., between 19 and 31 
nm for a single layer and 3d nn < £ < 4.2d nn (13 — 19 nm) 
for a multilayer. The upper and lower limits of £ corre- 
spond to, respectively, the lower and upper limits of the 
freezing temperatures T g that are compatible with the ac 
susceptibility experiments described above. Our magnetic 
memory experiments point to significantly shorter corre- 
lation lengths £ ~ 4.4 nm. 

Conclusions. — The central result of the present 
study is that we observe the same magnetic memory and 
critical slowing down in two-dimensional nanoparticle ar- 
rays, as well as in multilayers, suggesting that the un- 
derlying physical behavior is also the same. The slowing 
down of the ac (and dc) susceptibility curves measured on 
all these samples can, in fact, be accounted for using the 
activated law [Eq.([2j] that is typical of two-dimensional 
spin-glasses. These results suggest, therefore, that the 
glassy magnetic dynamics observed in these materials is 
associated with a phase transition occurring at T g = 0, 
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Fig. 6: (Color online). Zeeman energies estimated from the 
relaxation of ZFC magnetization curves measured after aging 
the sample at zero field and at T w = 15.8 K for a waiting time 
tw = 10 4 s, before the application of a magnetic field H. The 
inset shows the time dependence of the magnetic viscosity, de- 
fined as dM/d log(t) with M being the sample's magnetization, 
measured for (from bottom to top curves) H = 5 15, 20, 30, 
50, 60, and 80 Oe. This quantity shows maxima, marked by 
grey dots, at the effective age of the system that decreases with 
increasing H . The dotted line is a guide to the eye. 



rather than with a conventional spin-glass transition with 
a finite T g . This conclusion is supported by the results of 
magnetic memory experiments, which show that dynami- 
cal magnetic correlations are rather short ranged and, in 
any case, shorter than expected for canonical spin glasses. 

As mentioned in the introductory section above, the 
observed facts disagree with the prediction, derived from 
Monte Carlo simulations, that Ising-like spins interacting 
via dipole-dipole interactions should undergo a spin-glass 
transition below a finite temperature [17] • Establishing 
the origin for this discrepancy is beyond the scope of the 
present work. Here, we content ourselves with discussing 
possible deviations of real materials from the ideal con- 
ditions set by such models. In our opinion, an important 
aspect to be considered is the, unavoidable, distribution in 
particle sizes. In our multilayers, the distribution is prop- 
erly described by a Gaussian function of width ajj ~ 0.7 
nm, which is equivalent to roughly ±1 atomic layer and 
provides an indication of the good homogeneity of these 
samples. This narrow size distribution leads, however, to 
an extremely large dispersion in the relaxation times r 
associated with the magnetic anisotropy of the nanoparti- 
cles. Using the parameters U and To given above (see also 
the second reference in [TO] for further details), it follows 
that intrinsic timescales separated by more than 13 orders 
of magnitude can coexist at temperatures near or below 
Tb- A possible consequence of this enormous dispersion 
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is the following. Smaller, and therefore faster, relaxing 
nanoparticles are able to immediately react to spin flips of 
larger (slower) ones, minimizing their mutual interaction 
energy. We have previously shown that this effect accounts 
for the modification of the average relaxation times by in- 
teractions, at temperatures close to [T5]. We might 
speculate with the possibility that the disorder in relax- 
ation times also hinders the growth of magnetic correla- 
tions at lower temperatures. For instance, the formation 
of negatively polarized magnetic clouds surrounding the 
largest nanoparticles can screen dipolar interactions be- 
tween them. Clearly, further theoretical studies that in- 
clude effects of disorder and nonequilibrium dynamics are 
required to clarify the nature of the collective magnetic 
response in nanoparticle arrays. 
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